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1. INTRODUCTION 


The army is currently investigating methods for arming 
helicopters with weapons of greater firepower. A major diffi- 
culty stems from the relatively weak helicopter structure which 
cannot sustain the recoil of most large caliber weapons. Al- 
though recoillessrifles obviate one problem, the blast field 
created by propellant gases discharged through the rifle nozzle 
creates other problems. Damaging peak overpressures of several 
pounds per square inch would be produced on some paris of the 
helicopter fuselage. 


Constraints which arise from considerations of loading, 
firing, etc., prevent mounting the rifle so that the blast field 
does not affect the helicopter's structure. Moreover, severe 
weight penalties are involved in shielding the fuselage or 
employing ducts to carry away the nozzle blast gases. The net 

8 effect of these constraints is to make the feasibility of arming 
helicopters with recoilless rifles depend on finding a means for 
attenuating the intensity of the blast field. 
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A particular concept for attenuating the blast field was 
investigated. The concept involves attaching a short cylinder 
filled with water or solid particles to the rifle nozzle to 
delay the emergence of propellant gas. The flowing propellant 
gas would be slowed by driving the particles* out of the cylinder, 
thereby reducing the blast pressure, if the expelted liquid or 
solid does not give up its acquired momentum to the atmosphere 
too rapidly after it emerges from the cylinder. 


The approach taken in this investigation is to develop a 
computational method and program for describing the blast field 
produced by the particles expelled from the cylinder. This is 
accomplished by extending the numerical methods presently used 
for solving unsteady compressible single component flows to a 
multiphase flow where the velocity and rasperature of the particles 
can be different from that of the gas in which they are suspended. 
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BLAST ATTENUATION 
a. Basic Mechanism 


The blast field behind a recoilless rifle depends mainly 
upon two factors, the total energy of the propellant gas flowing 
through the nozzle and the duration of the flow. Since the 
comparative importance of these factors varies with the distance 
from the nozzle, it is desirable to have an assessment of the 
distances in which one or the other has a predominant influence. 
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The distribution 2f pressure, density, energy or velocity 
in the blast field can be regarded as a system of traveling waves 
in which the local wave speed is the speed of sound. As the pro- 
pellant gas flows into the blast field and raises the temperature 
by the effect of compression or by the convection of gases of 
higher temperature, thc vave speed increases. A consequence of 

this is that a profile of pressure versus time that formerly 

rose gradually to the peak pressure in the vicinity of the nozzle, 
changes so that it rises more rapidly at greater distances. As 
one considers increasing distances from the nozzle the pressure 
profile eventually steepens to fourm a shock wave. The process 
takes place as illustrated in Figs. la and lb. For distances 
sarge enough for the profile to steepen up to a shock wave the 
pressure is independent of the duration of the nozzle flow and 
only depends on the total energy of the pas and the distance from 
Jw nozzle. Lf the shock pressure were still strong compared to 
itis ambient pressure at these distances it would approach the 
Yaylor's result for a point source explosion (Ref. 1), 


po B,/R (1) 


de Taylor, G., "The Formation of a Blast Wave by a Very Intense 
Explosion,’ Vol. 201. 
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where 


E, = total energy released by explosion 
or in propellant gases 


p * pressure 
R = realistic distance from source . 


Although Taylor's result may not be useful for a quantitative 
estimate cf the blast pressure in the problem of concern, the 
indicated linear dependence of the pressure on the total gas 
energy is probably a good guide. 


These preliminary considerations of the propagation of the 
blast wave show that for positions in the blast field beyond 
the distance where a shock wave forms, the only means for 
attenuating blast intensity is to reduce the total energy in 
the propellant gas. At lesser distances, nozzle flow duration 
is a factor and prolonging it would result in reduced peak 
pressure. An approximation of the distance required for the 
formation of the shock wave can be made as follows. The tem- 
perature of the propellant gas increases the speed of sound 
two to three times over that which exists in ambient air. The 
peak chamber pressure is reached about 6 msec after ignition. 
Equating the travel distances for waves emitted at the initial 
and peak pressure conditions and assuming that the waves at the 
peak condition travel at a speed three times greater than the 
others, yields a distance of 10 ft. 


The area of interest on the helicopter extends from about 
3 to 10 ft from the fuselage. For much of this area it is 
apparent that both of the factors discussed will influence the 
blast pressure, but at the extremity the peak pressure is 
determined by the total energy in the propellant gases. 


b. Approaches for Implementing the Blast 
Attenuating Mechanism 


Two recoilless rifle configurations were considered in which 
particles are expelled. Each epitomizes a different blast attenua- 
tion mechanism: (1) reduction of total energy, and (2) proionga- 


tion of flow duration. 
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(1) Davis Gas 


A rifle based on this principle would consist of a straight 
tube in which particles are expelled from one end ard the pro- 
jectile from the other. Since the momentum of the particles must 
equal that of the projectile to fulfill the requirements for 
cancellation of recoil, the tube length required is large unless 
the weight of the particles is much larger than that of the pro- 
jectile. For particles of weight equal to the projectile, the 
tube length would be twice that of a conventional rifle and, 
therefore, impractical. ‘The blast pressure is reduced because 
the gas does not escape from the rifle freely until the particles 
are expelled from the tube. An estimate of the energy removed 
from the propellant gas by accelerating the projectile and 
particles may be obtained by assuming that the expansion is 


isentropic, 


e 
= = (volume ratio)” 1, 0.66 (2) 
1 
where 
e, * initial internal energy 
ey, = internal energy when propellant gas emerges 
v = ratio of specific heat, 1.21 for propellant 


gas volume ratio = 1/7, chamber-to-chamber 
plus barrel volume 


About one-third of the energy in the propellant gas is removed 
and the approximate effect on the blast field is proportional to 
the value given for the ratio of the emergent to the initial 


internal energy. 


This concept provides a means for reducing the blast intensity 


and shows that it is feasible to attenuate the blast pressure by 
expelling particles, The question remaining is, do the expelled 


particles release their acquired momentum to the air after exiting? 
If the particles remain as a coherent mass, they would not; if they 
expand rapidly, the drag forces acting on the particles would cause 
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them to transmit their momentum to the air again and produce a 
blast field. This problem is dealt with by mesns of the two- 
dimensional blast field code discussed in the following section. 
A Listing of one- and two-dimensional codes is provided in 
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Aprendices A and B. 


(2) Short Cylinder Concept 

This concept was developed in.an attempt to achieve blast 
attenuation without the excessive length of the Davis gun. The 
cylinder from which particles are expelled is shorter and a 
nozzle is used between the cylinder and rifle chamber to permit 
‘4u- chamber pressure to ‘build up properly. With the shorter 


cylinder, the amount of energy that can be removed from the 
propellant gas is less than with the Davis Gas concept. To have: 
a significant reduction: in the intensity of ‘the blast field, 

tle effect of flow duration has tobe exploited. In this study, 
the particles are not compacted and ‘the propellent gas is 
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permitted to leak through void spaces, thus prolonging the flow 
period. Since the particles are not attached. to the cylinder 

and do not transmit an axial load :to it, the. recoilless proper- 
ties of the rifle are not affected significantly by the cylinder... 


}. MATHEMATICAL MODELING 


The physical problem to be modeled involves the motion of 
liquid droplets or solid particles through a gas in which strong 
compressible effects take place that include the formation of 
shock waves. Because the temperature and velocity of the particles 
are not usually the same as of the gas at the same location, 
various interactions take place between the gas and the suspended 
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particles. The effects of drag and heat transfer cause inomentum 
and energy to be transferred between the gas and the particles. 
the motion and temperature of the gas therefore, are different 
than in an analogous flow case without particles. 


A numerical method is developed for the solution of this prob- 
lem that is patterned after the same techniques used for the 
numerical solution of single component flows. Conservation rela- 
tions for the mass, momentum and energy are developed for the gas 
and particle phases and these relations are then expressed in 
finite difference form for solution by numerical methods. Ap- 
proaches to the development of the conservation relations and some 
sulution for particular flow’ cases are summarized by Soo (Ref. 2) 
and Harble (Ref. 3). 

a. Conservation Relations 

The conservation relations are developed for an elemental 
cell as indicated in Fig. 2. Particle sizes and mean separation 
distances are assumed to be small relative to the size of the 
cell. ‘This limits the applicability of the method to problems in 
which the mean particle separation distance is smali relative to 
the system size and the scale of the phenomena of interest. 


(1) Mass Conservation 


Equating the accumulation of gas or particles to the net 
flux of gas or particles into the cell and the contributions due 
to evaporation and other effects, gives the following expressions 
for the gaseous and particle phases: 


2. Soo, S. L., Fluid Dynamics of Multiphase System, Blaudell 
Publishing Co., (196/). 


3. Marble, F. E., “Dynamics of Dusty Gases," Annual Rev. 
Fluid Mech., Vol. 2, (1970). 
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where 
Cy = concentration of gas 


cy = concentration of matter present as particles 


N = number of particles per unit volume 
t = time 
ty = velocity of particles 


b 


= velocity of gas 


= ¢c 
® a 


= mass loss of a single particle by vaporization 


W... = mass addition by propellant burning 


(2) Momentum Conservation 


When one considers the motion of a single particle the governing 
equation is obtained by equating the Lagrangian acceleration to the 
drag force and the pressure difference across the particle. 


ok 


u ~~ 
d as « - 
™ pe D- V Vp (6) 
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where | 
D = drag force vector 


3 
L] 


mass of the particle 


= pressure 


<= Jv 


= volume of the particle 
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For the purpose of the computational method to be developed, a 
Eulerian expression (i.e., at a fixed point for the derivative), 
is required and can be obtained if the number of particles within 
the cell indicated in Fig. 2 is large enough so that the mean 
velocity is reasonably well defined and can be assumed to be 


continuous. 
se a 
ms t= -mu,-V U, -D- Vv Vp (7) 


Multiplying through by the number of particles and using Eq. (4) 
gives the following for the mOmentum of the particle phase. 


BE 7 7 Vitg Uys Uy - ND - 5 Vp o 


where p, is the density of the material of which the particles 
are composed, 


Although the dyadic notation used (V-cu; u) is inconvenient, 
it permits the development of a conservative computational scheme, 
i.e., one in which the accumulation of momentum in the computational 
grid is exactly consistent with the fluxes of momentum over the 
grid boundaries. The alternative and more common expression for 
acceleration (u‘V u), does not permit this. 


The momentum equation for the gas phase may be obtained by 
subtracting the equation for the particle momentum derived here 
from the equations for total momentum of both the gas and the par- 
ticles. However, the experience obtained with various different 
computation procedures with this problem, showed that a more stable 
program is obtained if the computational procedure is based on 
the equations for the total momentum of both the gas and particles 
and the momentum of the particles alone rather than on the equations 
for the momentum of the gas and particle phases, ee 
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Equating the accumulation of the momentum to the fluxes across 
the cell boundaries and the effects of pressure, gives the following 
expression for the total momentum of both phases, 


ie (c, u, + Cy oy) we Ve(c 


g Ug? Ye te, us uy) - Vp (9) 


Terms relating to particle drag do not appear in this equation since 
they control only the transfer of momentum between the gas and 
particle phases and thus do not add to or substract from the momentum 


sum. 


(3) Energy Conservation 
Considering again a single particle, the following equation is 


obtained for the temperature by equating the rate of change of 
stored sensible heat to the losses by convection and evaporative 
cooling. 


DT, 
mM Cog pe * 7 (Q, + Q,) (10) 


where 


Cpa = specific heat of the particle 


Q. = cooling rate of particle by convection 


Q. = cooling rate of particle by evaporation 


An Eulerian representation fcr the rate of temperature change can 
be obtained when the number of particles within the cell is large 
enough to permit a reasonable definition of a mean particle 


temperature, 
de 4Ty _ 
Cod Tas + Vecy Uy Ty / fo. N(Q, + Q,) (11) 


An expression for the total energy in both phases is obtained by 
equating the rate of accumulation of energy to the fluxes and the 
work done by pressure effects at the cell boundaries. 
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where B 
e, = internal energy of the gas phase e 

VF = void fraction Fi 

8 


As in the case ‘of the momentum equation, the computation scheme 
is more stable if the computations are based on the total energy 
and the energy in the particle stream. The gas phase energy is. then 
obtained by calculating the difference in these quantities rather 
than by calculating the gas energy from an equation for the energy 
of the gas phase. 


b. Subsidiary Relations 


To obtain solutions for the above set of equations, it is 
necessary to connect the concentration, pressure and density with 
an equation of state, to eveluate the drag forces and particle heat 
transfer rates that transfer momentum and energy between the gas 
and particle phases. 4A propellant burning law is also required 
to predict the rate at which mass and energy is added to the 


propellant gases. : 
Equation of State 4 2 

Tabulated values for the thermodynamic properties of the | 4 

propellant gas from 23 different military propellants are given z 3 

in Ref. 4. For computational purposes § polynomial representation i 3 

was used to Fit the equation of state data. The following nine 4 4 

term polynomial was fitted by least squares to the data for M8 i 2 

= i propellant with a maximum discrepancy of only 0.3 of 1 percent: a 4 
2 = 
3 a 

4. Baer, P. G. and Bryson, K. R., Tables of Gomputed Thermo- “he 3 
denice. Properties of *Military un Propellants, BRL.Memo. “Fog 3 
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where 
= density 


Py 


e, = internal energy, including energy of 
& formation as defined for e in Ref. 4 


Numerical values for the coefficients A; are given in Appendix A 
in the subroutine entitled EQSTAT of the one-dimensional code, 


The gas density to be used in this equation of state is 
cetermined from the amount of gas existing in a unit volume and 
the void fraction. 


4 ; (14) 
p, =e /V (15) 


@ Drag Forces 


Dray forces acting on the drops are calculated From the drag 
coefficient and the dynamic pressure as follows: 


D=C, tae ju, -u,| (ay > 4.) (16) 
dG Py SAIN Seeenemenen 


where 
a #* particle diameter 


C; * drag coefficient 
p, * gas density 
Absolute values of the velocity difference are used in the manner 


shown te preserve the proper sign and direction of the drag force 
irrespective of the sign of the velocity difference. 


4. Baer, P. G. and Bryson, K. R., Tables of Computed Thermo- 
dynamics Properties of Military Cun Fropeliante. BRL Memo. 
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Coefficients of drag Cz are; obtained from the particle Rayne }es 
_ number Re as shown, Ref. 5 ; 


Cy = at _ Re «< 2.05 . , oF 


a 0.6 2.05 « Re « 486 


0.44 486% Re 


Woe re 


re ears 
pa fu, - 
ge ed | 
- : ug s * 
iy = gas viscosity 


_@ Heat Exchange 


whe rate at which heat is transferred from the gas to the particulate 
phase is obtained from correlations for the heat transfer coeificient oO 
on a single particle (Ref. 6), 


ha = 2 +0,34 Re’*® p_ 1/3 (203 
é : r : 


h = heat transfer coefficient 


x 
n 


thermal conductivity of the gas 


* 
a 


Prandtl number of the gas 


The heat transfer rate, Q, is given by 
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Q =a” h(t, = Ty) . (21) 


5. Perry, J. H., Chemical zi i yeer n Handb: Kk, McGraw-Hill 
Publishing Co. - RTSS0F. 


6. McAdams, Heat Transmission, McGraw-Hill Publishing Co. (1954). 


cd * 


ae the gas temperature. 


@ Propellant Burning Law 


The rate at which gas and.energy is added to the proper tant 
yas will be computed from surface regression rate of the bornin 
propellant. Watervliet Arsenal provided the following correlation 
tor the propellant to be considered. 


R = 0.00186 p?°83 


where p is the pressure in pounds per square inch and R is the 
regression rate of the propellant surface in inches per second. 


c. Coniputational Procedure 

For the purpose of presenting the numerical method by which 
solutions to the preceeding equation may be obtained, it is 
useful to use a more compact notation. Each of the differential 


Eqs. (3), (4), (5), (8), .(9), (11) and (12) can be expressed in 
the following form. 


orn ar(uq. + &) oy or(wt + e_) (23) 
n n n n n Mes oe 
Xe” 7 TUF wR er Ge Ten 
where r is the radial distance from the axis and z is the 
distance along the axis. Each of the variables a to g are 
defined in Table I. The velocity components, u and w, are of 
the phase considered. 


The numerical procedure is based on a finite difference form 
for each of the equations represented by Eq. 23. Conceptually, 
the field of interest is ¢’vided into descrete cells and the finite 
difference equations descrive the rate of accumulation of mass, 
number of particles, momentum, energy, etc. in terms of the fluxes 
of these quantities across the cell boundaries and the other 
quantities appearing in Eq. 23. In the particular vrid system 
used in this study, velocity is defined at the cell boundaries 
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and particle number, concentration, pressure energy and temperature 
are defined at the cell centers as indicated in Fig. 3. ‘The fluxes 
across the cell boundaries, are computed in accordance with the 

donor cell concept described by Gentry and Martin (Ref. 7). That is, 
the flux of particles for example, is calculated from the particle 
density in the cell which donates the particles. 


Because velocity, in the computational procedure used here, 
is defined at the cell boundaries,computations for the momentum 
of the gas and of the particles are based on special cells displaced 
one half the distance between the grid lines so that it is centered 
over the points where the velocity is defined. The momentum, pu, 
is thus based on the velocity at the special cell centers and an 
interpolated value for the density at that point. Similarly, the 
kinetic energy at the centerof the cells where internal energy is 
defined is based on interpolated values for the kinetic energy. 


Our experience with this computational procedure is that 
spurious negative values for the internal energy and pressure are 
obtained less frequently than with the more usual procedure in 
which the velocity is calculated for the same positions as the 
density and energy. 


For the study of the blast field behind the weapon, the method 
sould be implemented in a two-dimensional grid as shown in Fig. 4. 
The grid spacing increases with distance from the nozzle to accommo- 
date a large region without an excessive number of grid points. 


_Two computer programs were prepared for solving the blast field 


problem: (1) a one-dimensional program which describes the interior 
ballistics of the rifle and the flow from the cylinder; (2) a two- 
dimensional program which uses the flow from the cylinder ag input 
and calculates the blast field. 


7. Gentry, R. A.; Martin, R. E.; and Daly, B. J., "An Eulerian 
Differencing Method for Unsteady Compressible Flow Problems," 
Computational Phys., Vol. 1, pp 87-118 (1966). 
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Figure 3) ELEMENT SPATIAL DISTRIBUTION WITH INTERLEAVING 
Oi VELOCITY FOINTS BETWEEN POINTS AT WHICH 
DENSITY, ENERGY, PRESSURE AND COMPCSITION ARE 
CaALCUIA TED 
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4. RESULTS 


The interior ballistics and nozzle blast field were calcu- 
lated for the conditions defined by the following values: 


charge 8.1 1b 
web 0,061 in. 


I 
ent 


i iat 


i 


Ee bore 105 mm 
pe barrel length 140 in. 
EF: chamber volume 200 in. 
= cylinder diameter 7.4 in. 

cylinder length 24 in. 


throat diameter (effective) 3.26 in. 


Results for three different conditions in the attached cylinder are 


presented: 


e no water in the cylinder 
e 5 1b of water with a drop size of 5 mm 
@ 5 1b of water with a drop size of 1 mm 


Results obtained with the one-dimensional model which furnished 
the inputs for the two-dimensional model are given in tie following 
figures. Chamber pressures are given in Figs. 5-7. 


These were obtained by calculating the average value of the 
pressures at nodes in the chamber region of the gun. A value of | 
approximately 8000 1b was obtained for the peak chamber pressure 
in each case. The peak pressure is shown to increase slightly when 
water is used in the cylinder and when the drop size is diminished. 
This is due to the resistance to the flow or propellant gas developed 
by the water drops which diminish the leak rate of the rifle. The 
abrupt drop in chamber pressure that occurs at about 15 msec, 
coincides with the burnout of propellant. Values obtained for the 
projectile travel and velocity while in the barrel are shown in 


Figs . 8-10. 
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The output of the one-dimensional program was used as an 
input condition for the two-dimensional program. Early time 
examples of the spatial pressure distribution are shown in 
Figs. 11-12. The flow out of the attached cylinder is directed 
to the right and the cylinder exit is located on the horizontal 
axis in the center of the smallest contour. As one would expect, 


the pressure profile is elongated in the direction of flow. For 


the purpose of evaluating the effect on the blast pressure of water 
in the cylinder attached to the rifle, the pressure was listed for 
selected points for each case. The location of these points A-F 

is indicated in Fig. 13, and the pressure pulse profiles in 


Figs. 14-16. These pressure pulse profiles show that the peak 


pressure increases as water is added to the cylinder and the drop 
size is diminished. 
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5. SUMMARY OF RESULTS AND CONCLUSIONS 


Contrary to the reasons advanced for expecting blast attenua- 
tion, the results obtained show that the peak pressure is increased 
rather than reduced by expelling water. The duration of the pulse 
is shorter, however, so that the total blast load represented in 
the pulse may be somewhat reduced. The most probable explanation 
involves the effect of flow duration and perhaps also the ability 
of the particles to transfer momentum to the blast field after 
expulsion from the cylinder. . 


In the study of the short cylinder concept it was thought 
that an optimum particle size might exist that would prolong flow 
duration, thereby reducing the peak blast pressure. It appears, 
however, that particles of suitable size, on the order of 1 mn, 
also have a drag-to-weight ratio so large that the momentum absorbed 
from the propellant gas in the cylinder is given back to the air 
when the particles emerge. The remaining hope for reducing the 
blast pressure with a concept involving particle expulsion is to 
use considerably more water so as to absorb more energy from the 
propellant gas, or to exploit other phenomena not investigated in 
this study, such as quenching the gas by the cooling effect of 
droplet vaporization. 
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APPENDIX A 
ONE-DIMENSIONAL GAS PARTICLE FLOW 
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